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Abstract 

We investigate the evolution of the bound state of negatively charged massive 
particles (CHAMPs) with light elements and discuss its cosmological consequences 
and the constraint. By numerically solving the Boltzmann equation, we study the 
time evolutions of such bound states. Since most of negative CHAMPs are captured 
by 4 He, its bound state is positively charged and couples with the electromagnetic 
plasma. When charged particles constitute a dominant non-relativistic component, 
density fluctuations of matter cannot grow due to the acoustic damping. This results 
in the suppression of matter power spectrum from which a severe constraint can be 
obtained. By arguing constraints from other aspects of cosmology, we show that the 
constraint from large scale structure gives most stringent one in some representative 
cases. 



1 Introduction 



Long-lived charged massive particles (CHAMPs) can exist in various extensions of the 
standard model of particle physics such as supersymmetry (SUSY). One of such example is 
a slepton, a superpartner of leptons, which can be stable if it is the lightest supersymmetric 
particle (LSP) and R-parity is conserved. However, the abundance of such stable charged 
massive particles would be severely constrained [1], in particular, from experiments of the 
deep sea water [2-7]. Although there might have been a mechanism to prevent them being 
captured into the Earth [8] and in such a case the constraints may not be applicable, a 
scenario with stable CHAMPs would generally not be viable. However, in some other 
scenarios, CHAMPs are unstable, and they can constitute a dominant component of non- 
relativistic particles in the early Universe. For example, when the gravitino is the LSP, 
which can be easily realized in gauge-mediated SUSY breaking models [9-11], the next 
lightest supersymmetric particles (NLSP) may be CHAMPs, and they can be long-lived. 
Although such unstable CHAMPs can evade the constraint from the sea water, they affect 
other aspects of cosmology. It has been rigorously investigated that the decay of such 
massive particles would destroy light elements synthesized by big bang nucleosynthesis 
(BBN), from which we obtain constraints on the properties of CHAMPs such as the decay 
rate and its abundance [12-23]. Although the considerations of the decay also applies to 
neutral massive particles, there is another important effect on BBN which is specific to 
the charged particles: the bound-state effect. Negative CHAMPs can form a bound state 
with positively charged light elements, which affects the BBN reaction rates and their 
abundances [24-40]. In fact, as will be shown later, most negative CHAMPs (referred 
to as X~) are captured by 4 He, a double-positively charged element. If this bound state 
( 4 HeX~) is stable for some time in the course of the history of the universe, it would also 
affect other aspects of cosmology in addition to BBN such as large scale structure. Thus a 
detailed investigation of how the bound states are formed and evolve would be important 
and interesting, which is one of the main topics in this paper. 

When a particle possesses an electric charge before recombination, it is tightly cou- 
pled with the plasma of electrons and photons. Thus such charged particles, which are 
supposed to be the bound state ( 4 HeX~) here, can also participate in the acoustic oscilla- 
tions. When the bound state of negative CHAMPs constitutes a dominant component of 
non-relativistic particles, fluctuations of (non-relativistic) matter cannot grow due to the 
"acoustic damping" caused by the acoustic oscillations, which is different from the case 
with a standard neutral cold-dark matter (CDM) model. This results in the suppression of 
the matter power spectrum at some scales [41-43]. Thus the bound state ( 4 HeX~) should 
have a significant effect on large scale structure, and the consideration of this issue can 
place a unique bound on the properties of (negative) CHAMPs, which is another topic we 
are going to focus in this paper. 

The organization of this paper is as follows. In the next section, we first carefully 
investigate the evolution of the bound states of X~ with some light elements by numerically 
solving the Boltzmann equation. Then we discuss the effect of the (charged) bound state 
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on large scale structure and the damping of matter power spectrum. In Section HI we 
discuss the constraint on the property of CHAMPs from large scale structure and some 
other aspects such as BBN and the CMB spectrum. The final section is devoted to the 
summary of this paper. 

Unless otherwise stated, throughout this paper rii, rrii and Y$ denote the number den- 
sity, the mass and the yield variable (= rij/s with s the entropy density) of a particle "i", 
respectively. 

2 Evolution of the bound state 

We are interested in the bound-state formation of X~ with a light element, which occurs 
after the cosmic temperature becomes lower than 30 keV. Until that time, most of the 
standard BBN processes should have almost been finished. Under this circumstance, the 
Boltzmann equations for the time-evolution of the number density of bound states (NiX~), 
denoted as n^ NiX -), with iVj = p and 4 He are expressed by 



dt 



+ (o~bnd, 4 Het') 
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( m v m x T \ 3 ^ 2 _ F , T 
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-((r cx v)(n4 Hc - n { 4 HcX -))n {pX -), (2) 

where n x - is the number density of free X~, T x is the decay width of X, E b 4 Hc ~ 
337.33 keV and E bp ~ 24.97 keV are the binding energies of ( 4 HeX~) and (pX~) [31], 
and n-tHe and n p are the number densities of 4 He and proton including both free and 
bound-states. The masses of the bound states are given by m^ Ke x-) = m^ + mx — E b 4 Hc 
and m(px-) — m p + m x ~ E bp , respectively. The thermally-averaged recombination cross 
sections (cr bnd 4 He v) and (crbnd.pw) for the 4 He and p bound-states formation are given by [27] 

Knd, 4 Hc^) ~ 98.46 aE ^ He (3) 
m Rc Vm4 Ue T 

(a hndtP v) ~ 24.62 (4) 
where T is the cosmic temperature, and a is the fine structure constant. 



2 



The terms in the third lines of the right-hand side of Eqs. ([I]) and (j2J) represent a 
charge-exchange reaction, 



(pX~) + 4 He -> ( 4 HeX~) + p, (5) 

with its thermally-averaged cross section being denoted as (cr cx v). This process is effective 
just after (jpX~) has been formed. Recently it has been reported that a rate of this charge- 
exchange reaction is more rapid than the Hubble expansion rate [44] at the formation epoch 
of (pX~) which is approximately given by 



H 



where we took the yield variable of free 4 He to be V4 Hejfree = 5.4 x 10~ 12 . This means that 
the produced (pX~) is immediately destroyed, and rather X~ is included into ( 4 HeX~) 
after the destruction. 

In Fig. [U we plot the time evolution of /bnd,Ni = Y(N i x~)/YN i fl with Ni = 4 He (left) and 
Ni = p (right), where Yn u o and Y x - are the initial values of and Y x -, respectively. 
Here the initial value of Y^ t means the one after iVj's standard BBN processes have finished 
and well before its bound-state formation starts. Y x - is the initial value of Yx- well before 
X~ decays, and/or its bound-state is formed. The dashed lines denote solutions of Saha's 



equation^ 1 



Y iNiX - } = (^f^)~ i/2 ^ /T ^ ~ Y^XXx- - Y {NiX - } ). (7) 

From the figure, we can see that the bound states form at around T ~ 10 keV and 1 keV 
for ( 4 HeX~) and (pX~), respectively. 

The time-evolution of the ( 4 HeX~) formation excellently agrees with the solution of 
the Saha's equation. On the other hand, deviations from the Saha's equation can be seen 
in case of (pX~). The behavior of the time-evolution of (pX~) seen in Fig. [1] (right) can 
be easily understood as follows. If Y x - >0 3> Y4 He0 , sufficient amounts of free X~ exist, 
independently of the detail of ( 4 HeX~) formation, and the abundance of (pX~) approx- 
imately follows the solution of Saha's equation. On the other hand, if Y x - >0 l4 Hc0 , 
the abundance of (pX~) should deviate from Saha's equation. Furthermore, the charge- 
exchange reaction becomes important when ((Jbnd,pv)Y x -Y p ~ (cr cx v)Y( J) x~)Y4 He . Then the 
formation of (pX~) is balanced between the production and the destruction processes, and 
the abundance becomes approximately the order of Y( pX -) ~ 0.5 x 10~ 5 Yx-,o(^c/keV) -1 / 2 . 
Thus the production of (pX~) stops at around T c ~ 0.6 keV for Yx- o < 10 -12 seen in 



* x Note that Saha's equation represents the equilibrium between X~ and iVj and assumes only one 
component of Ni. Therefore, there should be deviations from this solution when we consider two compo- 
nents of Ni(= p and 4 He). However because the dashed lines still give us useful information to understand 
behaviors of the numerical solutions, we also plot them in Figs. Q] and [2] 
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Figure 1: Time evolution of /bnd,Ni = Y(N l x-)/YN l fi with Ni = 4 He (left) and Ni = p 
(right), where Y^o and Y x - are the initial values of and Y x -, respectively. The 
dashed lines denote solutions of Saha's equation. In these figures we assumed that X~ is 
stable. 



Fig. [T] (right). This feature is consistent with Fig. 2.7 in Ref. [45] and what was stated in 
Ref. [44]. The reason why X~s are included mainly into 4 He is that the number density 
of X~ is much smaller than that of electron, wihch is completely different from the case 
of the standard recombination of electron. 

On the other hand, in Fig. [2] we plot the time evolution of fx,N l = Yr NiX -\/Y x - j0 with 
Ni = 4 He (left) and N = p (right). Similarly to FigJT]the dashed lines show the result of 
Saha's equation. From Fig. [2] we see that most of X~ are captured into the bound state 
with 4 He if Y x - o is smaller than 10~ 12 . This means that (pX~) disappears immediately 
after its formation for Y x - q ^5 10~ 12 . Then we find that almost all X~s are captured by 
4 He, and form the bound state ( 4 HeX~). Because the abundances of other singly-charged 
nuclei such as deuterium and tritium, are much smaller than that of proton, we can omit 
contributions from deuterium and tritium. This result tells us that the bound states of 
X~ cannot be neutralized even after {pX~) could have been formed. 

The fact that the total electric charges of the dominant non-relativistic components 
cannot be shielded has a strong impact on the structure formation. We discuss this issue 
in the next section. 



3 Effect on Large scale structure 

In this section, we discuss the effect of the charged bound state on large scale structure. 
In fact, the discussion below also applies to charged (massive) particles themselves when 
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Figure 2: Time evolution of fx,N, = Y{NiX-)/Yx-,o with Ni = 4 He (left) and Ni = p 
(right). The meanings of the labels and lines are the same as those in Fig{TJ 



they constitute a dominant non-relativistic component. Thus in the following, we use 
"CHAMPs" to indicate both the free CHAMPs and the bound states of CHAMPs with 
light elements which have a net electric charge. When there exist long-lived CHAMPs, 
they couple with photon-baryon fluid, thus matter density fluctuations oscillate under 
the scales which enters the horizon before CHAMPs decay. Thus fluctuations of matter 
cannot grow, then the matter power spectrum is suppressed on corresponding scales, which 
is called "acoustic damping" in literature. 

An explicit calculation has been done on how the matter power spectrum is suppressed 
in [41,42]. To obtain the constraint on the decay rate rigorously, we need to compare the 
matter power spectrum with observational data. However, when CHAMPs are dominant 
component of non-relativistic matter, the matter spectrum is abruptly suppressed at the 
damping scale. Thus the evaluation of the damping scale would be enough to obtain 
the constraint on CHAMPs. Hence in the following, we simply make an estimate of the 
acoustic damping scale as a function of the decay rate. 

The scale under which matter power spectrum is suppressed, which we denote kx 
in the following, can be estimated as follows. For fluctuations of the scale which enters 
the horizon before CHAMPs decay (t < tx where Tx is the lifetime of CHAMPs), they 
are damped by the acoustic oscillation. On the other hand, fluctuations of the scale 
which enters the horizon when t > Tx, in other words, they can be assumed to be usual 
neutral dark matter, fluctuations of such scales grow with time as in the standard casa^ 2 . 

# 2 Here we do not consider free-streaming of nonthcrmally-produced dark-matter particles because it 
depends on the kinetic energy of the dark matter just after its production by the X~ decay, which would 
be highly model-dependent. To obtain a conservative constraint, wc neglect this effect here. 
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Since the horizon-crossing occurs when k = aH with H being the Hubble parameter, the 
characteristic scale kx under which the matter power spectrum is suppressed is defined by 

k x = aH\ t=Tx . (8) 

We assume that the CHAMPs decay during radiation-dominated (RD) epoch. Then H 
is related to the cosmic time as H — l/2t. On the other hand, during RD, H can be 
also written as H 2 ~ p r ad/3Mpj = H^VL^d'^ . Putting these together, we obtain the 
characteristic scale kx as 

To make some rough estimate, we recall that fi ra( j/i 2 ~ 4.15 x 10~ 5 and H ~ h/(3 x 10 17 ) s, 
we obtain 

k x ^ 10%/^Mpc" 1 . (10) 
V r x 

For example, when Tx ~ 1 s, the corresponding damping scale is k x l ~ 0.1 kpc. For 
Tx ~ 10 6 s, the damping scale becomes k x x ~ 10 2 kpc, which is the order of the galaxy 
scale. Here it should be mentioned that the damping scale down to k x l ~ 1 kpc can 
be probed with future observations of QSO-galaxy strong lens system [43]. Thus the 
constraints from large scale structure would be much more stringent and such future 
observations give us a lot of information on CHAMPs. 



4 Constraints on CHAMPs 

Now we discuss the constraint on the abundance and the decay rate of CHAMPs from some 
cosmological observations, paying particular attention to that from large scale structure. 
We show the constraint on the Y x - vs. Tx plane, where Y x - is the yield variable of 
X~ which is the number density of X~ to the entropy density ratio. In the following, 
we assume that the energy density of CHAMPs is fixed to give the present dark matter 
density if they are stable, (i.e., the energy density of CHAMPs before the decay is fixed 
by this requirement.) Thus the yield variable and the mass are related as 

y x - * 4 x io- 12 n DM h 2 (1^1). (n) 

V rn x - J 

In Fig. [31 we draw the acoustic damping scales k x x = 1 Mpc, 100 kpc, 10 kpc, 1 kpc 
and 0.1 kpc by the vertical long-dashed lines. If we confirm structures larger than a 
scale, we can exclude the right region of the corresponding line. Here it should be noted 
that unstable CHAMPs can erase the structure of small scales, and thus can solve some 
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Figure 3: Constraints on the yield variable and decay rate of negative CHAMPs. The 
energy density of CHAMPs is fixed to give the present dark matter density if they are 
stable and then the mass is also fixed by this requirement. The corresponding masses 
are given by the horizontal dotted lines with the label near the lines (See also Eq. (TTTT) ) . 
The right side of vertical long-dashed line is excluded by the requirements that 0.1 kpc, 
1 kpc, 10 kpc, 100 kpc, and 1 Mpc, size structure should not be erased, which is plotted 
from left to right. 100 kpc corresponds to typical galaxy size structure. Upper region of 
thick solid line is excluded by 6 Li overproduction by the Catalyzed BBN of the bound- 
state effect [26,30-34]. Notice that the constraint from CBBN does not depend on the 
CHAMP mass. We also show the constraint from BBN and CMB spectrum with the 
mass of CHAMPs being fixed. These constraints should be interpreted with some care 
in this figure since the mass is fixed independently of the requirement of Eq. ( TTTT) for 
these cases. Upper region of thin solid line is excluded by the BBN constraints from 
an minimal assumption of possible photodissociation for the visible energy of the decay, 
E vis = 10 3 GeV and 10 6 GeV, respectively. They are obtained by appropriately scaling 
the result of [17,20]. Upper region of dashed lines are excluded by fi- and y-distortions 
of CMB spectrum for same visible energies. Their labels are located near the lines. For 
the case where the lifetime is longer than the age of the Universe, the constraint from the 
deep sea water may apply, which is shown as the shaded region. 
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problems regarding the discrepancies between observations [46, 47] and predictions of N 
body simulation in ACDM model [48-52] 



#3 

Another important constraint comes from the overproduction of e Li by the Catalyzed 
BBN through ( 4 HeX") + D -> 6 Li + X~ induced by the formation of ( 4 HeX~) [26]. 
Observational fraction of 6 Li/ Li gives a upper bound on Yx- [26,30-34]. Notice that the 
constraint from CBBN does not depend on the CHAMP mass. The upper bound on Yx- 
for Tx > 10 6 sec is Y x - < 2 x 10~ 15 . Thus, requiring the abundance of CHAMPs gives the 
present dark matter density if they are stable, the constraint from large scale structure 
becomes more relevant than that from CBBN for nix ^ 10 6 GeV. 

For reference, we also show other possible constraints placed on Yx- and tx from 
the photodissociation of BBN [17] and CMB y and /x-distortion of CMB spectrum [56] 
with the mass of CHAMPs being fixed independent of the relation of Eq. (11 II) . Since 
the way of fixing the mass is different from that in other constraints, these constraints 
should be interpreted with some care in the figure. Here we assumed that X~ emits 
electromagnetic particles and take some representative values for the energy injected from 
the decay of X~ as E vis = 10 3 GeV and 10 6 GeV. The lines of the constraints from the 
photodissociation and the CMB are obtained by appropriately scaling the result of [17,20] 
and [56], respectively. In addition, here we have assumed that the branching ratio into 
electromagnetic particles (£> v is = E v - ls /mx) would be the order of unity. The readers 
can easily obtain these bounds from photodissociation and CMB distortions by scaling 
the branching ratio into electromagnetic particles correspondingly. When the lifetime of 
CHAMPs is longer than the age of the Universe, the constraint from the sea water may 
apply [2-7]. In Fig. [3j we fixed the energy density of CHAMPs to become the same as 
the present DM density. Thus the constraint from the sea water is Yx- < 4 x 10~ 18 , 
which corresponds to m>x = 10 8 GeV. As mentioned above, the constraint from large scale 
structure does not depends on the mass and the yield variable. From the figure, we can 
see that large scale structure can give the stringent constraint in some parameter regions. 

So far we have discussed the case that the dominant component of the nonrelativistic 
particles in the universe could only be the negative CHAMPs. It would be trivial that 
long-lived positively-charged particles with tx ~ 10 6_8 s can be simply excluded by the 
same reason from the viewpoint of the large-scale structure without taking into account 
neither the bound-state formation nor their neutralization if they become the dominant 
component. 



# 3 It has been discussed that non-thermal production of warm dark matter can also erase the small 
scale structure by the long free-streaming length due to its relatively large velocity dispersion. [43,53] 

Furthermore, see also Refs. [54,55] and references therein for another idea to erase small scale structures 
by introducing bound-state formation through hidden gauge interactions of dark matter. Note that in 
their models the formation of the bound-state means the kinetic decoupling and the end of the acoustic 
oscillation of the hidden-charged dark matter. If they considered hidden 4 Hc additionally as well as 
standard cosmology, their situations might be changed. 
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5 Summary 



In this paper, we have investigated the evolution of the bound state of CHAMPs with 
light elements and have shown that the negatively-charged massive particle cannot be 
neutralized even after its bound-state formation with protons at around 0.5 keV. This is 
because the charge-exchange reaction by free 4 He through (pX~) + 4 He — > ( 4 HeA~) + p 
is much more rapid than the cosmic-expansion rate, and almost all X~ will be included 
into ( 4 HeX _ ) for Yx-,o~ 1CT 12 , which is positively-charged. 

This gives a high impact on the formation of the large-scale structure if those charged 
particles are dominant non-relativistic components of the universe like cold dark matter 
at the cosmic time 10 6 s. Then any galaxies cannot be formed by the suppression of 
the density perturbation through the acoustic oscillations. This simply means that the 
lifetime of the negative CHAMPs should be tx < 10 6 s at longest. 

As was discussed in the text, future observations of QSO-galaxy strong lens system can 
probe the structure down to k x l ~ 1 kpc. Those future observations will reveal the nature 
of the long-lived CHAMPs. In this letter, we have presented the result for the case where 
the abundance of CHAMPs is fixed to be the present-day dark matter density if they are 
stable. However, it would be interesting to investigate the case of the energy density of 
CHAMPs being changed, which will be the issue of a separate paper. 
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